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I. INTRODUCTION 

Monolayer and multilayer fine-particle deposits at sur- 
faces are of interest for a wide range of applications, in- 
cluding photonic crystals, quantum dots, heterogeneous 
catalysts, sensors, and microarrays 0, H, H, 0, IS S 0- 
There has been a recent drive to explore and utilize 
particles smaller than the traditional colloid-size (few 
microns to sub-micron), namely down to nanoparticles 
(means dimensions of 0.01 of a micron and smaller, i.e., 
sizes of order 10 nm and below). Quantification of the 
kinetics of synthesis, aggregation, and surface interac- 
tions of nanoparticles requires new experimental probes, 
but also new theoretical techniques. Furthermore, the 
surfaces with which fine particles interact, can now be 
pre-patterned to control and modify the particle at- 
tachment kinetics and the resulting deposit morphol- 
ogy 0, S H ©, El, EL E El- Presently, experiments 
have produced patterns on the submicron scale, but the 
trend is down to nanoscale. Thus, deposition kinetics will 
be on an artificially formed lattice or another pattern, 
which can improve catalytic and reactivity properties of 
the particle-covered final surface. From a theoretical per- 
spective, such processes pose interesting challenges, in- 
cluding identification of the parameters that control the 
properties of the resulting deposit. 

In the present work, we report a detailed study of the 
influence of a pattern consisting of square shaped cells 
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in which centers of circular particles can land (e.g., pro- 
jections of spherical fine particles depositing in a mono- 
layer). The cells are in turn arranged on a planar sub- 
strate in a square lattice array. We consider the pro- 
cess of random sequential, fully irreversible adsorption 
of fixed size particles (disks). As usual for random se- 
quential adsorption (RSA) processes [li . HH Il6l li~7l . [iH 
El US El, EL EL HI HI Hi, we assume that particles 
cannot overlap: the arriving disks, randomly and uni- 
formly transported to the surface, are deposited only if 
they do not overlap earlier deposited ones (and provided 
their centers fall within the square pattern of the allowed- 
deposition cells). Deposition attempts of disks that do 
not satisfy these conditions are rejected, and the parti- 
cles are assumed discarded (transported away from the 
surface). 

Thus the RSA model [U, El, EL El, El, EI as- 
sumes that the effects of the particle-particle interactions 
and particle-substrate interactions can be accounted for, 
approximately, by purely geometrical restrictions and 
features. The excluded volume constraint represents 
particle-particle interaction which is assumed short-range 
repulsion on length scales shorter than the particle sizes. 
Particle-substrate interaction is assumed to result in ir- 
reversible binding on the time scale of the process. Fur- 
thermore, the details of the particle transport to (and, for 
rejected particles, away from) the surface are lumped into 
the assumption of uniform flux of deposition attempts 
per unit surface area. In studies of RSA, one is inter- 
ested in characterizing the jammed state morphology at 
large times, when a dense RSA deposit is formed and 
no available particle landing sites are left, as well as the 
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approach to the jammed-state coverage. In this work we 
focus on the former aspect of the process: the jammed 
state properties. Despite its intrinsic simplicity, the RSA 
model provides a surprisingly rich set of limi ting behav- 
iors and morphologies 
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The paper is organized as follows: In Section |TlJ we 
start by defining the model, while in Section IIII1 we 
analyze, analytically and numerically, the jammed state 
properties for the case when the excluded-volumc inter- 
action is functional only within individual landing cells. 
For geometries for which the excluded- volume interaction 
extends beyond individual landing cells, extensive Monte 
Carlo simulations are reported in Section |IV] and |V] We 
argue that the jamming coverage is not sufficient to study 
the jammed state, particularly in this case. The radial 
particle-particle correlation function is considered as the 
property of interest, in Section |V1 which also offers some 
concluding remarks. 



II. THE MODEL 

Our primary goal is to study irreversible monolayer de- 
position of identical hard-core spherical particles on flat 
patterned substrates. This is obviously equivalent to de- 
position of disks with the hard-core, "no overlap" exclu- 
sion on a plane. The particle centers are only allowed to 
adsorb within well-defined bounded regions. These land- 
ing cells will be for simplicity modelled as square "land- 
ing areas" allowed for the disks' centers. We point out 
that other cell shapes [1, E3, E3, 0, E E| are fea- 
sible and relevant, e.g., rectangular, circular, etc. More- 
over, we assume that the cells themselves are regularly 
positioned in a square lattice array, and, again, we note 
that other lattices are possible, e.g., the triangular lat- 
tice. The lattice structure has its own crystalographic 
unit cells, which are larger than the landing cells. 

In the present work, we consider particles of fixed ra- 
dius, tq. This is obviously a theoretical idealization. Ex- 
perimentally, the particle sizes and shapes will always 
have some dispersion. For dispersions above 0(10%) 
of the mean size, particles are considered polydispersed. 
However, syntheses of uniform spherical colloid and nano- 
size particles with polydispersity as low as 4% has been 
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62j . so that the monodisperse approximation 
is quite realistic for many systems of interest. 

Thus, we assume that identical particles arrive with 
flux F, i.e., that the rate, per unit time, of deposition 
attempts of disk centers at the substrate is F per unit 
area. The particles arrive to the surface randomly and 
uniformly. A deposition attempt fails if the disk's center 
falls outside the allowed-landing cells, or if the arriving 
disk overlaps one or more previously adsorbed ones, the 
latter mimicking the excluded volume interaction. These 
are the assumptions typical of the RSA model, and the 
binding is assumed irreversible: the particles do not de- 
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FIG. 1: Dashed lines delineate the square lattice unit cells. 
Adsorption of disks can only take place when their centers fall 
inside the solid squares, i.e., the allowed-landing cells. 



tach from or move (diffuse) on the substrate on time 
scales of the dense deposit formation. This latter as- 
sumption is usually a very good approximation for col- 
loid deposition [24[,[38||, but can be questioned, e.g., for 
protein deposition: indeed, studies of RSA-type models 
with particle rearran gem ent on the surface have also been 
reported [H 011, LzH ■ In the irre- 
versible RSA model, the deposit density initially grows 
linearly with time, t. However, as the particle density 
increases, the hard-core exclusion leads to slow-down of 
the deposition process. Ultimately, for large times the 
"jammed state" is approached, at a density lower than 
that of close packing and with no long-range order in the 
particle positioning, but with no gaps left for any addi- 
tional particle deposition. 

For purposes of our modeling, we assume square 
landing-cells of linear size a, positioned in a square lat- 
tice with unit-cell size a + b (obviously, a, b > 0). This is, 
of course, an idealization, but we find that this geometry 
already offers a rich pattern of deposit morphologies. We 
consider RSA of disks of radius r . Without loss of gen- 
erality, we rescale the substrate lengths relative to the 
diameter of the disks. Specifically, we define 
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We comment that in studies of RSA the time is also usu- 
ally rescaled, by a factor inversely proportional to the 
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FIG. 2: A disk fails adsorption onto the substrate because (a) 
its center does not fall within an allowed-landing cell, or (b) 
it overlaps with a previously adsorbed disk. 



rate and particle "volume," here (irr^F)" 1 . Figure [T] il- 
lustrates the geometry of the lattice (in dimcnsionlcss 
units). 

The deposition rules are illustrated in Fig. [2] Specifi- 
cally, we do not allow particles to attach on top of each 
other, so that multilayer deposits are not formed [72], [73| . 
Multilayer particle deposition was studied in various ex- 
perimental and theoretical contexts [H, [71, IzS III, . 
In this regard our present model represents a gen eralized 
version of lattice RSA El, M S3, S El H IzS 
HO, l8l|. o n par with such gen eralizations as RSA of mix- 
tures [Til, m m m mm m mm m m 13 or 

deposition on finite-size substrates 14 1. 



III. THE JAMMED STATE 

In the present model, the "phase diagram" is in the 
space of the two parameters, (a,/J), the size of the cells 
and the distance between neighboring cells, in reduced 
units. The most obvious quantity to consider is the cov- 
erage, measured by the fraction of the total surface cov- 
ered by particles, 0(t), in the jammed state for various 
values of (a,/3). The jamming coverage, 9j, is a prop- 
erty of the state in which no further particles can adsorb, 
obtained in the large-time, t, limit, 



h = lim 6{t) 

t — ► CO 



(3) 



If the landing cells are too far apart from their nearest 
neighbors, then particles at different cells will not be able 
to "see" each other through excluded volume. Specifi- 
cally, for cells at a distance (3 > 1 apart from each other, 
particles (disks) attempting adsorption cannot overlap 
other disks in different cells. We denote this as the non- 
interacting cell-cell adsorption (NICCA); see Fig. [3] For 
distances between cells j3 < 1 a particle attempting ad- 
sorption can overlap with a previously adsorbed one be- 
longing to a different cell, thus leading to a failed deposi- 
tion attempt. We denote this as the interacting cell-cell 
adsorption (ICCA). 

One can choose the cell size to have a maximum, pre- 
determined close-packed number of particles inside it. 
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FIG. 3: The major subdivisions in the two-parameter space. 
For cell-cell separation (3 < 1 we have the interacting cell- 
cell adsorption (ICCA), while for j3 > 1 we have the non- 
interacting cell-cell adsorption (NICCA). For cell sizes a < 
1/V% we have a single-particle-per-cell adsorption (SPCA), 
while for a > 1/V2 we have multiparticle-per-cell adsorption 
(MPCA). Limiting cases of the model, as well as subregions 
with interesting properties, are discussed in the text. 



For a < at most a single particle can adsorb inside 

any given cell. We denote this as single-particle-per-cell 
adsorption (SPCA); see Fig [3] For cells with a > 1/V2, 
more than a single particle can fit in the cell, and we 
denote this as multiparticle-per-cell adsorption (MPCA). 

For the remainder of this section, we consider the 
NICCA case defined by (3 > 1, which implies that a disk 
attempting adsorption with its center landing in a par- 
ticular cell will never overlap a previously adsorbed disk 
in another cell. Thus, the global kinetics of deposition 
decouples into independent local kinetics at each land- 
ing cell. Therefore, for this range of (3 values the model 
is equivalent to continuum RSA on finite-size substrates, 
with somewhat unusual boundary conditions that parti- 
cles (disks) can "stick out" of the finite a x a square as 
long as their centers are within the square. The morphol- 
ogy and other physical quantities of interest of the global 
jammed state are determined by the jammed states for 
finite system sizes for the perscribed a value. 

Let us consider the a and (3 values to have NICCA 
with SPCA, where the latter case holds for a < 1/V2; 
see Fig. U for a typical configuration. The kinetics corre- 
sponds to that of lattice RSA of monomers, since each cell 
is certain to have a single particle at the jammed state. 
The difference relative to the lattice RSA, is in the par- 
ticle positions, which here are uncertain within the order 
of the size of the cell. While perhaps theoretically least 
interesting, such patterning provides for the most "con- 
trolled" particle deposition in applications. Since each 
cell ends up having a single particle, the jammed-state 
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FIG. 4: (Color online) Typical configuration of a region of 
25 x 25 unit cells, for a = 0.6 and f3 = 1.2, in the jammed 
state. This snapshot corresponds to the NICCA-SPCA upper- 
left region in Fig. [3] 



coverage is simply 



FIG. 5: (Color online) Configuration of a region of 25 x 25 unit 
cells, for a = 1.2 and (3 = 1.2, in the jammed state. Particles 
attempting adsorption never overlap previously adsorbed ones 
in different cells, but contrary to Fig. [4] each cell can now 
adsorb more than one particle. Since j3 = 1.2 > 1, the kinetics 
of adsorption at each cell is decoupled from that at other cells. 



4(a + /3) 2 



(4) 



which holds for p > 1 and < a < l/y/2. Finally, 
for a = and j3 > 1, we recover the true lattice RSA 
of monomers which are of unit-radius disks deposited at 
the centers of square cells (3 X j3. 

In the regime of NICCA with MPCA, a possibility 
opens up for having two or more adsorbed particles in 
each cell; see Fig. [5] for a typical configuration. Since 
the kinetics of deposition is decoupled in the sense ex- 
plained above, one can make some conclusions regarding 
the coverages for a given number of particles. 

Let us consider the case of up-to-two particles per 
cell as an example. It is obtained for values of a in 
the range 1/-/2 < a < (1 + v / 3)/2v / 2- The coverage 
for the close-packed (maximally packed) configurations 
is given by n/2(a + (3) 2 for two particles per cell and for 
P > 1. The closed-packed coverage will change discon- 
tinuously at a = (1 + v / 3)/2v / 2, as illustrated in Fig. [5] 
However, our numerical results, also shown in Fig. [6l in- 
dicate that the value of the jamming coverage for the 
RSA model, which for more than one particle per cell is 
less than the maximal coverage, remains continuous at 
a = (1 + V3)/2V2 and also at a = l/y/2. This behav- 
ior continues for larger number of particles per cell, see 
Fig. [5] In the close-packed problem, the highest coverage, 
at fixed /?, occurs for the smallest cell, i.e., the smallest 



TABLE I: The table shows the jammed coverage values in the 
NICCA-MPCA case, i.e., for /? > 1 and a > l/y/2. In the 
first column, n is the number of particles per unit cell of the 
close-packed situation. The values of a n , defined in the text, 
are shown in the second column. The third column gives the 
close-packed coverage values for the range of up-to-n particles 
per cell. In the fourth column, the jamming coverage values, 
6 j, from simulations are presented. These values are for a 
representative choice /3 = /?* = 1.2 (see text), and for the 
a values shown in the second column. Finally, in the fifth 
column, the error estimate, aj, in the jamming coverages is 
given. 
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value of a = a n allowing the prescribed number of par- 
ticles, n. Figure [6] illustrates the discontinuous coverage 
increments at a n , followed by a decrease oc (a + ft)~ 2 ■ 

The simulated results for 9j, for a < l/y/2 follow the 
close-packed values, since the two problems coincide in 
this a range. However, at a = the RSA coverage 
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FIG. 6: (Color online) Coverage values obtained by Monte 
Carlo simulation of our RSA model (solid line, the jammed- 
state coverage) and by direct calculation for the close-packed 
configurations (dashed line), both for (3 = 1. Notice the dis- 
continuities in the values of the close-packed coverage as op- 
posed to the smooth variation in the RSA case. The jammed- 
state coverages for /3 = 2, calculated according to relation 
are also shown for comparison (dotted line). 

is continuous, since the probability of having a second 
adsorbed particle in any given cell remains small for a 
values above, but close to a-i = 1/V2. This property 
reflects the stochastic nature of the RSA model, i.e., in 
most cases the first particle adsorbs at a position inside 
the cell that blocks the chance for the second particle to 
adsorb later. A similar observation applies as a crosses 

Q!n=3,4,...- 

Our numerical results for 9j were actually obtained for 
a representative j3 value, (3* = 1.2. The coverage for any 
other j3 > 1, for a given fixed value of a, can then be 
calculated from 

ej(a,^ = (^pj6j(a,n- (5) 
The error estimate for the coverage is given by 



TABLE II: Tabulation of numerically simulated RSA 
jammed-coverage values, at fixed /3 = /3* = 1.2, for the range 
of q from 0.72 to 1.40. The error estimates are also shown. 
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where the index i numbers the Monte Carlo runs, and N 
stands for the total number of runs. From this definition, 
after some algebra one can show that 
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The above observation and the appropriate parameter 
and coverage values are summarized in Tables HI which 
in particular gives jamming coverage values 6j(a n , /?*). 



The simulated jamming coverage values for varying 
a are given in Table |TTJ In the present RSA prob- 
lem, for j3 = (3* = 1.2 the minimum of the cover- 
age is 6 JmiB = 0.21199 ± 0.00006, and it occurs for 
a min = 0.745 ± 0.005. Finally, in the limit a — > oo one 
recovers the well-known RSA of disks on continuum sub- 
strates [n h H3, m, m n m, m m, m • our value 

for the coverage is 0.54706 ± 0.00008, which should be 
compared to the range 0.54700 ± 0.00006 recently esti- 
mated in [g^l. Note that the system size for which this 
new estimate of the jamming coverage was obtained was 
2048 x 2048 particle diameters. The respective close- 
packed value is n/2y/3, corresponding to the n = oo case 
in Table [J 
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IV. MONTE CARLO SIMULATIONS FOR 
INTERACTING CELL-CELL ADSORPTION 

In this section, we begin our consideration of the ICCA 
regime. In this case, correlation can develop beyond sin- 
gle cells, and therefore extensive Monte Carlo simula- 
tions were warranted. Here we consider the jammed state 
coverage and morphology snapshots. However, particle- 
particle correlations have to be considered for a fuller 
quantitative description. This is taken up in the next 
section. Counting system sizes as integer multiples of 
the unit cell side, a + (3, we took 500 x 500 as the finite 
system used in numerical simulations, unless otherwise 
stated. Periodic boundary conditions were applied both 
horizontally and vertically. 

A full description of the algorithm used, which repre- 
sents an extension of the one introduced in [U HH , will 
be detailed in a separate publication [93| . Here, we sum- 
marize relevant details that might be pertinent to the 
direct interpretation of the results. We used the adap- 
tive mesh-cell scheme [13, HH to increase algorithmic effi- 
ciency, particularly at late times. Our adaptive mesh-cell 
scheme allowed for mesh-cell subdivision to half of the 
size, when the number of available mesh-cells dropped be- 
low 75% of its original value. Such mesh-cell subdivisions 
allow to identify and exclude sub-cells already blocked by 
earlier deposited particles, thus improving the calculation 
speed of adsorption in the remaining "open" cells. One 
can show that this scheme does not bias the adsorption 
of particles (93|. We also comment that this approach 
can follow the time-dependent kinetics, though we only 
consider the jammed state in the present work. We also 
devised improved overlap tests of up to four surrounding 
(previously adsorbed) particles, which improved the algo- 
rithm performance. Another extension of the published 
algorithm [34l . [9l| , involved a natural implementation of 
the presence of the pattern, as detailed in [93|. 

The use of the algorithm allowed a detailed study of 
the structure of the jammed state. The number of re- 
alizations per simulation was 100. We report the values 
of coverage and the corresponding error estimates in Ta- 
bles HI HH and IIIIl for the various cases studied in the 
present work. 

Contrary to NICCA, in ICCA particles attempting ad- 
sorption at a particular cell, can overlap a particle previ- 
ously adsorbed in a different cell. The ICCA regime is ob- 
tained for values of /3 < 1. The coverage values obtained 
in our simulations are given in Table IIIIl As opposed to 
the NICCA case, in ICCA the kinetics of adsorption is 
no longer "decoupled," and particles, or clumps of parti- 
cles that belong to the same landing cell, not only follow 
the square positioning pattern of the landing cells but 
can also become correlated with particle in other cells. 
The resulting deposit morphology and degree of ordering 
will depend on the geometrical parameters, as well as on 
the fully irreversible nature of the adsorption of the RSA 
model. In fact, depending on the values of a and (3, parti- 
cles attempting deposition can overlap others that belong 



TABLE III: The table shows simulated values of the coverage, 
6 j, and the corresponding error estimates aj, for several val- 
ues of a and (3, in the ICCA-SPCA and ICCA-MPCA regions 
of the "phase diagram" shown in Fig. [3] 



a 


P 


9j 


a B] x 10 4 


0.1 


0.14 


0.5377 


13 


0.2 


0.2 


0.50866 


7.6 


0.2 


0.28 


0.54865 


6.7 


0.2 


0.5 


0.46563 


5.0 


0.2 


0.8 


0.54499 


3.1 


0.3 


0.42 


0.48338 


5.2 


0.4 


0.2 


52428 


5.2 


0.4 


0.56 


0.54519 


3.9 


0.5 


0.7 


0.54326 


0.6 


0.6 


0.2 


0.53444 


3.9 


0.8 


0.2 


0.55103 


2.9 


0.8 


0.8 


0.33358 


1.5 


1.0 


0.2 


0.53529 


2.5 


1.0 


0.8 


0.39362 


2.4 


1.2 


0.2 


0.53725 


2.3 


1.2 


0.8 


0.45497 


2.4 



to cells more distant than the nearest neighbor cells of 
the landing cell. Therefore, in addition to the jammed- 
state coverage, a more detailed study in the ICCA regime 
should also involve consideration of particle-particle cor- 
relations in the jammed state, as addressed in the next 
section. 

Combining ICCA with SPCA, which holds for param- 
eter values (3 < 1 and a < l/v2 and corresponds to the 
lower left region in Fig. [3l can lead to particle configu- 
rations correlated beyond the cell pattern. We identify 
some limiting cases of interest. One such limit occurs for 
(3 < 1 and a — > 0, which correspond to a well-defined 
square lattice structure for the deposition of disk cen- 
ters. Since each adsorbed particle effectively "shades" a 
circle of unit radius, which is larger than the lattice con- 
stant (3 < 1, deposition does not correspond to that of 
monomers as compared to the NICCA-SPCA case dis- 
cussed above, as the adsorbed particle will surely block 
neighboring cells, and possibly more remote ones, de- 
pending on the value of (3. However, taking (3 = 0, the 
system no longer has a pattern, regardless of the value of 
a. Thus, the corner of the phase diagram near the point 
(a, (3) = (0, 0) is special, but we did not consider this re- 
gion because it is more mathematically interesting than 
physically relevant: see recent literature on the kinetics 
of this type dHH. 

Due to the interaction between particles at different 
cells, during the deposition process correlations can de- 
velop, resulting in nontrivial local particle arrangements. 
These correlations stem from the excluded volume inter- 
action between particles upon adsorption and they can 
induce local semi-crystalline order. Visually, see Fig. 
the "crystallites" are oriented along the diagonal direc- 
tion of the square lattice of square cells. However, it is 
well known [23,[!l|,[!2i that the RSA process alone cannot 
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FIG. 7: (Color online) Configuration of a region of 25 x 25 unit 
cells, for q = 0.2 and (3 = 0.5, in the jammed state. A parti- 
cle attempting adsorption can overlap a previously adsorbed 
particle in a different cell. This excluded volume interaction 
is responsible for correlations which result in locally-diagonal, 
semi-ordered domains as seen in this snapshot. 



FIG. 8: (Color online) Configuration of a region of 25 x 25 
unit cells, for a = 1.2 and j3 — 0.2, in the jammed state. 
For this low j3 value, the probability of a particle attempt- 
ing adsorption to overlap with one in a neighboring cell is 
appreciable, thus building up a somewhat longer-range diag- 
onal semi-ordering than that seen for the parameter values of 

pig. m 



impose long-range ordering. Indeed, the particle corre- 
lations in RSA are usually rather short-range. We find 
that the patterning of the surface can induce ordering 
over several lattice spacings, which reflects the symmetry 
of the underlying pattern. More generally, the ordering 
should also depend on the shape of the deposited objects 
[33l . l74j . Thus, the pattern does influence the creation 
of ordered structures in an otherwise uniform deposition 
process. However, the stochastic RSA dynamics tends to 
prevent the long-range order of the pattern from being 
fully "imprinted" in the deposited particle configuration, 
as observed in Fig. [7J 

Combining ICCA with MPCA, with the parameter val- 
ues (3 < 1 and a > 1/ v2, corresponds to the lower right 
region of the "phase diagram" in Fig. [3] In this regime, 
the excluded volume interaction can lead to deposit mor- 
phology with semi-ordering beyond a single cell. How- 
ever, the overlap with particles in neighboring cells can 
extend at most up to the second-nearest-neighbor cells 
(diagonally neighboring cells). In fact, for 1 > > l/v2 
(with a > l/\/2), the overlap can be at most up to 
the nearest-neighbor cells. As a increases, one expects 
a lower impact of the cell-cell excluded volume interac- 
tion on the morphology. Consequently, for more cell-cell 
exclusion effects, we need smaller cell sizes, with limited 
number of particles (but at least two) allowed in each. 
The cell-cell exclusion leads to a further reduction of the 
average cell population, is illustrated in Fig. [5J where 
a = 1.2 and (3 = 0.2. For these values of the parameters, 



each cell have enough area to accommodate up to four 
particles, but excluded volume interaction due to nearest- 
neighbor cells for this low value of (3 substantially lowers 
the average cell population, as compared, e.g., with that 
of Fig. [5] for a — 1.2 and f3 = 1.2: the average number 
of particles per cell in the case of a = 1.2 and j3 = 1.2 
is 2.487 ± 0.001, while that of a = 1.2 and (3 = 0.2 is 
1.3407 ± 0.0006. However, the case of = 0.2 still has 
a higher value of the coverage, 0.5373 ± 0.0002, while for 
/3 = 1.2 the value of the coverage is lower, 0.3391±0.0001, 
because of more void space between the cells. 



V. INTERPARTICLE DISTRIBUTION 
FUNCTION 

To further characterize the jammed state, we introduce 
the distribution function of the distances, r, between the 
centers of the adsorbed particles, P ra A (a, (3; r). The sys- 
tem is translationally invariant in terms of the integer 
multiples of the unit cell size (a + (3) , because we use pe- 
riodic boundary conditions both horizontally and verti- 
cally. Therefore, only displacement vectors between par- 
ticle centers matter when studying particle-particle cor- 
relations. For convenience, in order to avoid discussion of 
the "connected" part vs. the full correlation function, we 
normalized the correlations by counting only distances 
between pairs of particles within a cutoff distance, R. 



8 



The distances r and R will be assumed dimensionless, 
measured in units of the particle diameter. We found 
it appropriate to limit our study to separations up to 
R = 5(a+/3). Since our results were not accurate enough 
to study possible weak singularities that can develop in 



This distribution function is normalized as follows, 

f P rad (a,/3;r)rdr = l. (9) 
Jo 

The shape of the distribution function in the jammed 
state depends on the values of a and f3. The position 
of the first peak measures typical distances between the 
closest particles. To better understand the role of a 
and (3, we considered three families of distribution func- 
tions, as defined in the following subsections. The time- 
dependent kinetics of a distribution function will be stud- 
ied in a forthcoming article [93[. 



A. Effect of varying a on the distribution function 

We start by studying the influence of varying a on the 
semi-ordering of the jammed state. We carried out a sc- 
ries of simulations at fixed (3 = 1.2 and varied a. The 
results are shown in Fig. [9] The dimensionless center- 
center distance, i the horizontal axis of the plot — 

was further rescaled in terms of the unit cell size, to 
r/(a + 0), see Fig. [5] Keeping (3 > 1, here (3 = 1.2, 
corresponds to the NICCA case. In the NICCA-SPCA 
case, the first peak appears at a unit (rescaled) distance, 
since the distance to the closest particle, on average, is 
that to the nearest-neighbor cell. For NICCA-SPCA, 
well-defined peaks also appear that correspond to other 
underlying lattice distances defined by the square-lattice 
pattern. This is apparent in the distribution functions for 
a = 0.2, 0.4, shown in Fig. Ufa), with peaks at distances 
of 1.0, \/2 « 1.4, 2.0, and \fb « 2.2, etc. Increasing the 
value of a in the NICCA-SPCA regime increases the un- 
certainty in the position of the particle within the cell, 
i.e., it leads to peak broadening. 

Now, in the NICCA-MPCA case, the position of the 
first one or more peaks will depend on the value of a for 
a given value of (3 = 1.2 > 1; also shown in Fig. HJa). 
Additional peaks will reflect the intra-cell particle dis- 
tribution and can be positioned well below the unit-cell 
size. Peak broadening and peak-peak overlap are super- 
imposed in this case, but the pattern-induced tendency 



the jammed-state RSA correlations at particle-particle 
contact [111 [U, H2| , at r = 1 , and we were interested in 
the tendency for semi-ordering on length scales of several 
unit cells, we found it convenient to define 



(8) 

I 

for semi-ordering is still quite visible in the appropriate 
curves for a = 0.8, 1.4 in Fig. [UJa). 

In the ICCA case, the peak structure and the position 
of the first peak will be more complex due to the possibil- 
ity of overlap of an incoming particle with others belong- 
ing to neighboring cells. In terms of the distribution func- 
tion, the presence of the excluded volume interaction has 
an interesting effect of smoothing out the lattice-induced 
tendency for ordering, as shown in Fig. [9£b). The effect 
is particularly noticeable when one compares these plots 
with those for the same values of a in Fig.[H{a). The fea- 
tures of the distribution function seem to be determined 
primarily by the particle-particle jamming effects, rather 
than by the underlying landing-cell pattern. 

Finally, for f3 = 0.8, though the interaction between 
particles at different cells is present, it is not as promi- 
nent as for [3 = 0.2, and one observes an intermediate 
behavior, as shown in Fig. [H^c). The distribution func- 
tion is still smoothed out due to the jamming effects be- 
tween particles at different cells, but traces of the lattice- 
induced ordering remain, especially for a = 0.2,0.4. 



B. Effect of varying (3 on the distribution function 

In order to discuss the effect of the value of (3 on the 
structure of the jammed state, let us first consider fixed 
a = 0.2, with varying (3 = 0.2, 0.5, and 1.2, as shown 
in Fig. nWa). One observes that, as (3 increases, the 
distribution function becomes more detailed with peaks 
becoming sharper. There is also peak splitting, related 
to a lesser degree of excluded volume interaction between 
a particle attempting adsorption and another one from 
a different cell. We comment that for values of j3 > 1 
the general shape of the radial distribution function is 
no longer changing, since particles cannot overlap. 

The excluded volume interaction with particles belong- 
ing to neighboring cells also reduces the number of par- 
ticles effectively adsorbed in a cell as observed in snap- 
shots of the jammed state in Figs. [5] and El with a fixed 
at 1.2, while (3 changing from 1.2 to 0.2, respectively. 
For a — 1.2 and for increasing values of (3 in the interval 
(0, 1), one observes, see Fig. [TUJb), that the radial dis- 
tribution function reflects more structure from particle 



_ , r, \ Number of pairs of particles with distances in (r, r + dr) 

r A\CL 13' T) = — 

raa ' ' r dr (Total number of pairs of particles at distances < R) 
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1 1.0 



B= 0.20 
B= 0.50 
B= 1.20 
B= 1 .40 




a= 0.20, B= 0.28 
a= 0.30, B= 0.42 
a= 0.50, B= 0.70 
a= 0.60, B= 0.84 




FIG. 9: Plots of the distribution function for various values 
of a, at fixed values of /3: (a) /3 = 1.2, with a = 0.2, 0.4, 0.8, 
and 1.4; (b) (3 = 0.2, with a = 0.2, 0.4, and 1.2; (c) (3 = 0.8, 
with a = 0.2, 0.4, 0.6, and 1.2. 



FIG. 10: The distribution function (a) for constant a = 0.2, 
with P = 0.2, 0.5, 1.2; (b) for constant a = 1.2, with /3 = 0.2, 
0.4, 0.6, 0.8, 1.2; (c) along the diagonal line (0, 0)— (l/y/2, 1) 
in the ICCA-SPCA region of the (a, (3) "phase diagram," Fig. 
H with (a, /3) = (0.2,0.28), (0.3,0.42), (0.5,0.70), (0.6,0.84). 



C. Jammed-state structure along (3 = V^a 



arrangements inside the cells. As (3 increases, the posi- 
tion of the first peak also shifts to lower values. However, 
as far as the pattern-induced ordering is concerned, very 
little trace is left of it, and the curves are relatively flat, 
dominated by particle-particle jamming effects. 



As our last example, we consider the effect of chang- 
ing both a and (3 along the diagonal of the ICCA-SPCA 
region, which corresponds to values of < a < \j\f2 
and < j3 < 1. Here particles are larger than the cells 
and they can overlap particles from neighboring cells, 
not necessarily the nearest-neighbor ones. In the ICCA- 
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SPCA regime, particle deposition leads to highly corre- 
lated jammed structures. Though not specifically taken 
along the diagonal line, a snapshot of such a highly- 
correlated jammed state can be seen in Fig. As a 
rule, the tails of the distribution functions, see Fig.fTOFc). 
mostly coincide regardless of the values of both a and 
j3. This simply reflects the fact that these correlations 
are dominated by the excluded- volume jamming effects, 
rather than by the landing-cell pattern. However, corre- 
lations up to (scaled) distances of ~ 2.5 units do show 
parameter-dependent features. Specifically, the position 
of the first peak shifts to lower values of the distance, 
since the relative size of the particles compared to that 
of the unit cells decreases for increasing values of a (or 
0). 

D. Concluding remarks 

We studied RSA of disk-shaped particles on patterned 
substrates, specifically, square landing-cells positioned in 
a square lattice array, with focus on the jammed state 
properties. An efficient numerical algorithm was imple- 
mented to simulate the two-dimensional disk deposition. 
The interplay of the two geometrical parameters, the cell 
size and cell-cell separation, was found to have a striking 
influence on the deposit morphology and density, as well 
as on the in-cell particle population. We found that the 
deposit morphologies could be lattice like, locally homo- 
geneous, and locally ordered. The distribution function 
describes the degree to which the cell pattern affects the 
deposit morphology for various values of the cell size and 



cell-cell separation parameters. This ordering effect com- 
petes with the buildup of more random correlations due 
to the excluded volume interaction. 

Due to the increased use of particles of sizes well be- 
low 1 jtxm, future applications of deposition on patterned 
substrates will involve features at the nanoscale, which 
are quite difficult to manufacture with geometrical pre- 
cision and uniformity. Therefore, possible directions for 
future work should include consideration of the effects of 
randomness in the pattern, of particle diffusion within 
the cells, of particle detachment, as well as of differences 
in the interactions of the arriving particles with surface 
features, such as cell interior vs. edges. 
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